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Abstract 

The progress in nanofabrication, measurement technology, and meso- 
scopic transport theory has been expanding the held of shot noise. 
Although a wave-packet approach to DC shot noise of independent 
electrons at finite temperature was offered as an intuitive alternative 
to the sophisticated theories, actual shot noise data often behave more 
complicated than the derived simple expression. For example, so-called 
effective charge can deviate from elementary electronic charge due to 
correlated tunnelings. Also, there are cases where one wishes to know 
the full spectrum of the shot noise. It will be of great use if a handy 
method for the shot noise in various experimental situations is avail- 
able. In this article, a classical wave-packet approach to the shot noise 
is presented. The classical formulation provides a rigorous yet straight- 
forward formalism to compute the full spectrum and, furthermore, clar- 
ifies the structure of Fano factor and effective charge. Additionally, the 
role of realistic detectors and an application to cross correlation mea- 
surements are also discussed. The present method can serve as an 
intuitive complement to the full quantum mechanical and field theo- 
retical approaches. 
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1 Introduction 



1.1 Power spectrum of signal 

Suppose a signal I(t) is of our interest. When the signal spans a long time, 
the total energy of the signal, \I(t)\ 2 dt, may diverge. Hence, the power 
of the signal is more senseful. Imagine taking a segment of the signal I(t) 
and name it as I*r(t). 
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Here, let us imagine being a long enough time compared to any time scale 
of the system. Its Fourier transform is I,y{f)- Since I(t) is a real signal, 
I*5r(f) = I$"{~f)- We will often use lo = 2irf as well. We now consider the 
auto correlation function of I(t), C(t) = (I(t')I(t' + t)). 

,9^c 



(I(t')I(t' + t)> 



lim — / dt'lAt')lAt' + t) 

J J-&/2 
7 J- 3 



lim 



-3r/2 



Z7T Z7T 



1 



// 



lim 

2vr 2vr j7->oo ST 

dco du' n . 1 
hm — 

2vr 2vr S^oo S? 



du duj' . 1 
hm — 

2vr 2vr .sr^oo & 



sr/2 

e i(w'+w)t' 
i(uj + uj') 



dt > e i(u>+u)t> i^i^jy^'t 

9/2 



f\„ioj't 



,9/2 

1 
2 



UJ+LU' 

2 



dfdf' lim — 



■7 



iAu)i.Au'y Jt 



*(f + f')& 



IAf)lAf')e 2 * ift 



SI 



dfdf lim — 3?sinc(iT(f + f')£T) i<rtf)l<rU')j Jt 

s> -^too ,y L J 

dfdf hm ■!($(/ + f)lAf)lAf)e 2 * if,t 

3?—>-co ,y 



-1 



df lim l^(/)^(-/)e 2 -/* 



d/ lim 



l J ^(/)| 2 „27Ti/t 



(2) 



Hindoo is called the power spectral density S(f). The auto corre- 

lation function and the power spectral density are directly connected by the 
Fourier transform. This relation is called Wiener-Kintchin theorem. Since 
I(t) is a real function, both S(f) = S*(—f) contribute to the spectral den- 
sity at the frequency /. Hence, it is natural to redefine the shot noise power 
spectrum as S(f) = S(f) + S(-f). 

5,/) = ton Mi? (3) 

8>— >co iy 

1.2 Definition of shot noise 

The definition of the shot noise is the power spectrum of the current fluc- 
tuation, 51 (t) = I(t) — (I), appearing in the presence of driven currents. 
This is equivalent to the Fourier transform of the second order correlator of 
current fluctuation. The unit of the shot noise is A 2 /Hz. 

C(t) = (6I{t')6I{t' + t)) (4) 
S(f) = 2*[C(t)] (/>0) (5) 

In reality, the noise consists of several components such as external, in- 
strumentation, thermal, and shot noises. Especially, the shot noise is the 
component which appears in the presence of driven current. In the present 
article, we neglect the external, instrumentation noise. Therefore, the cur- 
rent fluctuation is nothing but the shot noise. Even at finite temperatures, 
as long as one focuses on the DC voltage bias range beyond the thermal 
energy, the formulation is directly applicable. Note that SI(t) is present 
even without external bias at an equilibrium finite temperature. This is 
the origin of the thermal noise, 4/cb0G, where G is the temperature and G 
is the conductance of the system, which is an example of the fluctuation- 
dissipation theorem. Here, the fluctuation is 5I(t) and the dissipation is 
G^ 1 . The thermal noise is known to exhibit a white power spectrum with a 
high-frequency cutoff at the order of 2fc ^ e 



1.3 What this method is about 

Martin and Landauer employed a wave-packet approach to calculate DC shot 
noise of finite temperature one-dimensional electrons impinged on a tunnel 
barrier [1 ] [2] [3j . Roughly speaking, the shot noise as a function of DC bias 
exhibits a V-shape with the slope proportional to the charge but rounded 
only in the small bias range comparable to the temperature. From this, one 
can extract the temperature and the charge of tunneling entity. Despite 
of its simplicity, it has been effective to explain various shot noise data 
even of highly correlated systems such as edge modes of fractional quantum 
Hall effects [3j[5][6]. Certainly, there are other sophisticated formalisms 
beyond Landauer's approach. Many experimentalists who just dived into 
the subject of electronic shot noise may find themselves a bit lost in front of 
literatures employing a heavy machinery such as Non-equilibrium Green's 
functions [7] [8]. Although it is a powerful method, computing the shot noise 
in various experimental conditions is often beyond experimentalists' tools 
(or impossible even for theorists). Furthermore, actual shot noise data can 
behave in more complicated manners than the simple V-shape. Sometimes 
one also wishes to know about the full noise spectrum. Hence, it is useful 
to have a handy method for the shot noise. 

This article provides probably the most elementary way to calculate 
the shot noise in various situations by employing a classical wave packet 
approach to the shot noise. The mathematical derivations is in a similar form 
as in the theory of signal processing [9] . The full expression of the shot noise 
contains two important terms, the power spectrum of the wave packet and 
Fano factor. The concept of effective charge will be also discussed. Despite 
of the simple formulation, it clarifies the structure of the power spectrum, 
Fano factor and effective charge in a trasparent manner. Therefore, this can 
serve as a complemental method to the full quantum mechanical or field 
theoretical formalism. 
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Figure 1: (a) An ordered train of wave packet, (b) Schematics of an exper- 
iment. The splitter may be quantum point contact, quantum dot, interfer- 
ometer, molecular bridge, etc. 

For clarity, we imagine a one-dimensional conductor and a splitter (parti- 
tioner /tunnel barrier) in the middle (e.g. quantum point contacts, quantum 
dots, etc). Current in such a system can be viewed as a train of wave packets 
(Fig. 1(a)). In an ordered train of wave packets, neighboring wave packets 
are separated by a time r. Each slot is either (multiply) occupied or empty. 
Important to note that the splitter does not fracture them. Each wave 
packet maintains its shape in the course of the propagation. Importantly, 
this method only deals with occupation of each slot and not the probability 
amplitude of the wave packets. As shown in the Fig. 1(b), the experiment 
consists of three step: (1) an ordered train of wave packets is driven from a 
source to the system; (2) the wave packets are partitioned in stochastic or 
correlated manners into two ways (e.g. transmitted and reflected, or right 
and left, depending on the realization); (3) each is deteced at detectors 1 and 
2 at the end of the system. Here, we start with an ideal detector, which has 
an infinite bandwidth and is able to resolve the shape of each wave packet 
and to check the occupation of each time slot. The case of realistic detectors 
will be also discussed. 
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2 Shot Noise Power Spectrum 

2.1 General form 

2.1.1 Current 

First of all, we need to describe current as a train of wave packets. Suppose 
that the wave packet shape is g(t), whose Fourier transform is JP[g(t)] = 
9(f)- 9(t) should be normalized as J g(t)dt = 1, which one can also regard 
as the quantum probability. The train can be expressed by the convolution 
of train of spikes and g(t). The time difference between the spikes is r, 
which is at the scale of r ~ -y. The fcth time slot appears at the time fcr. 
The expression of the train of spikes is as follows. 

N 

hi(t) = £ (t k = kr) (6) 

k=\ 

There are N slots. To describe the current after the system, we introduce 
the occupation of each slot, o& (k = 1 ■ ■ ■ N). 

N 

h(t) = J2o k S(t-t k ) (7) 
k=i 

Now the train of the wave packets is 

N 

h(t)*g(t) = (8) 
k=\ 

The electric current I(t) has to be multiplied with charge q. 

N 

I{t) = Y.l^-tk) (9) 

k=l 

This is absolutely everything one wishes to know about an electric current. 

2.1.2 Average current and current fluctuation 

Defining 2? = Nt, the average of the impinged current is 

N f *r/2 

< J *> = i im / 9(t-t k )dt 
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The shot noise is the specral density of the current fluctuation. The fluctu- 
ation in the occupations is given by 5o k = o k — {o k ). 

N 

5I(t) = Y,q6o k g(t-t k ) (12) 
fc=i 

2.1.3 Correlation function 



Now, let us calculate the correlation function. 
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(<5o ra+m Jo ra )^v_| m | is the average of So n+m So n over iV — \m\ events. The 
approximation from the 7th to the 8th line was done with one assumption. 
Namely, there are finite correlations, if at all, only among small m (m « 
N). Certainly, it is senseless that a tunneling event affects 1000 events later. 
Hence, {5o n+m 5o n ) N _\ m \ = (8o n+m 8o n ) in the large N limit. Therefore, the 
sum over m implicitly means summing over some small numbers (say, up to 
m « 10). We denote it as YTm without specifying the upper bound for m. 
Besides, we have defined the correlation function of g(t) as G(t). 

G(t) = J dt'g{t')g{t' + t) , G(f) = \g(f)\ 2 (14) 
2.1.4 Power spectrum 

The shot noise power spectrum is S(f) = S(f) + S(—f). Combining / and 
— / makes it enough to take into account only m > component. 



S(f) = q(Ii)^(So n+m 6o n )(e- imUT + e imuT )G(f) 

m>0 
* 

= 2q(Ii)^2(5o n+m 5o n )cos(mujT)G(f) 

m>0 

= 2q{o n )(I i ) {6 ° n + m * 0n) cos(muT)G(f) 

= 2q{o n )F(f)G(f){I i ) (15) 

Thus, the shot noise power spectrum is proportional to the charge q and the 
average transmitted current {I) and the frequency profile is embedded in 
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G(f) and F(f). The shape of the wave packets only comes in G(f) whose 
bandwidth is the order of t _1 . The Fano factor contains the information of 
the correlation among the tunneling events. 

2.1.5 Absence of white noise in ordered currents 

When 0^ = 1, the current is an ordered train of wave packets with no fluc- 
tuatioin, namely 5o^ = 0. Such ordered current is a signal with periodicity 
of r. Thus, power spectrum of such current has a disrete spectrum being 
finite only at the harmonics of r _1 . Hence, the white shot noise vanishes. In 
general, if Ok changes in any periodic manner with a periodicity p, its power 
spectrum is finite only at the harmonics of p" 1 . Any periodic (ordered) 
currents has no shot noise. 

2.2 Wave packet shape 
2.2.1 Case of spikes 

The simplest example of wave packet is spikes, delta function g(t) = S(t — 
tk). The Fourier transform, autocorrelation function, power spectrum, and 
resultant shot noise are 

&[8{t-t k )\ = 1 (16) 
G(t) = S(t) (17) 
G(f) = 1 (18) 
S(f) = 2qF(f)(I) (19) 

For the train of spikes, the power spectrum is white noise (suppose F(f) = 
constant) up to infinite frequency and only the Fano factor gives a frequency 
dependence. Of course, the infinitely sharp wave packet is not physical 
because one needs infinitely high energy to localize a wave into infinitely 
sharp. Regardless these concerns, it still captures the essence. 
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2.2.2 Case of square wave packets 

We consider the case of square wave packets of duration r\ and height 1/rj. 
g(t) = i for \t\ < \ and for \t\ > 2). The Fourier transform, autocorrela- 
tion function, power spectrum, and shot noise are 

*\g{t)\ = (20) 



1 



t 



V 
f_ 



G(t) = -l-U. (\t\<r,), (\t\> V ) (21) 



V 

G( f ) = r sm( : /7?) i 2 (22) 



S(f) = 2 q F(f)\ S -^^l] 2 (I) (23) 
L irfr] J 

For general r], the shot noise exhibits a white spectrum. Here, note that 
the case of rj = r means a strictly DC current withouth any gap between 
the wave packets. Regardig the section about the absence of shot noise 
in periodic signals, one can also understand it because the sine function 
becomes zero at every harmonics of frequency -. 



2.2.3 Case of sine wave packet 



We consider the case of sine wave packets g(t) = / ri ff . The Fourier 
transform, autocorrelation function, power spectrum, and shot noise are 

= 1 (l/l<y), (|/|>|) (24) 
G(t) = Sm( " /ot) (25) 

G(f) = 1 (|/|<f), (|/|>|) (26) 
S(f) = 2qF(f)G(f)(I) (27) 

In the limit of large /o, this is reduced to the case of spikes. Though it was 
presented here, it is not very physical actually. Since there is no negative 
probability in quantum wave packets. (The sine function takes negative 
value too.) In many transport situations, a DC bias eV is applied to the 
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system above the Fermi energy. One can view the electrons within the energy 
eV form wave packets (in the small enough bais range, one can assume the 
dispersion of electrons is linear ui = vpk, where vf is the Fermi velocity). 
h/eV is the time scale of inter-electron time difference. Hence, hvpjeV is 
more or less the electron's wave packet size. The power spetrum of such 
a current should have a high-frequency cut at the frequency eV/h as the 
highest energy of the system. The sine wave packets generates a rectangle 
power spectrum whose cut off is at / = ^. This gives an insight to the 
actual shape of the electron wave packets. 

2.2.4 Case of Gaussian wave packet 

- * 2 

We consider the case of a Gaussian wave packet g(t) = ^ 27TU e ' 
Fourier transform, autocorrelation function, power spectrum, and shot noise 
are 

&W)\ = y^-^ (28) 



G(f) = e- (wu)2 (30) 
S(f) = 2qF(f)e-^ 2 (I) (31) 



Gaussian wave packets are easy to handle and often employed. However, 
note that it has finite, yet vanishing, spectrum even beyond eV/h. 

We have seen several examples of the wave packet shapes. Each shape 
generates a unique power spectrum. However, there is one common feature. 
As far as the integral of the wave packet over time is unity (J dtg{t) = 
1), G(f) goes to unity in the limit of low frequency, G(f) — > 1. This is 
trivial because integrating over time is nothing but taking only the DC 
component. In the next section, we will see one example that G(f) lacks 
the low- frequency component. 
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2.2.5 Case of paired wave packets 



Next, consider two wave packets g'(t) which are temporally separated by rj 
and have opposite sign of charge, q and — q, whose sum is zero. Thus, this 
is an uncoventional situation, where finite shot noise without net current 
exists. 



g(t) = g'(t)- g'{t + V ) 
The autocorrelation of g(t) and its fourier spectrum is 

G(t) = [ dt' [g'(t) - g'(t + r,)] [g'(t + t') - g'(t + t' + r?) 

dt' [g'(t')g'(t + f) + g'(t' + V )g'(t + t' + V ) 

-g'ityit + t' + ^-g'it' + ^g'it + t') 
= 2G'(t)-G'{t + r])-G'(t-ri) 
G(f) = [2- e -»'-eHG'(/) = 4sin 2 (ojr))G'(f) 



(32) 



(33) 

(34) 
(35) 



Though the net current is zero ((/«) = (I) =0), we can define the particle 
(pair) number current such that 

= Jjm f = \ (36) 
(o k )N 1 (o fc ) 



(I P ) 



lim 
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r 



(37) 



,^->oo ^ r r 

Using this, shot noise power spectrum of random train of the pair is as the 

following. 

S(f) = 8g 2 (/ p )F(/) S m 2 M)G / (/) (38) 

The low-frequency component of the shot noise is vanishingly small. Al- 
though this example may sound unphysical, there is a model system that 
behaves as this, which is a capacitively-coupled ID conductors. Generating a 
wideband power in a conductor capacitively induces a high-frequency power 
with a low-frequency cutoff depending on the resistances and the capaci- 
tance of the counductors in the other conductor. The resultant spectrum of 
the other conductor resembles the power spectrum given above. 
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2.3 Fano factor and effective charge 

2.3.1 Case of Poissonian wave packets 

So far, we have not discussed o n in detail. The ordered train of wave packets 
is partitioned at the splitter and becoems a stochastic train of them. Quan- 
tities such as (o n ) and (o 2 ) can be calculated once a specific distribution is 
given. Let us consider the simplest example, Poissonian distribution. For a 
Poissonian process, the mean and variance are the same. 

{5oD = (o n ) F(f) = 1 (39) 

S(f) = 2q(o n )(I i )G(f) (40) 
5(/«0) = 2q(I) (41) 

Here, we used the fact that G(f ~ 0) ~ 1 at low frequencies / << r _1 and 
(I) = (o n )(Ii). The Schottkey formula for shot noise is obtained. 

2.3.2 Case of independent Fermionic wave packets 

Here, we assume that those wave packets follow the Fermi statistics. The 
impinged current is ordered such that o n = 1. The tunneling process is the 
binomial distribution. The wave packets are transmitted (reflected) with a 
probability T (1 — T = R) at each tunneling event. We assume there is no 
correlation among the events. 

N 

(N-n)\nl 

, where n is the number of transmitted wave packets. Given the distribution, 
the expectation values of n and n 2 are 

(n) = T^[T + R] N = NT (43) 
(n 2 ) = T—T—[T + R] N = (NT) 2 + NT(1 - T) (44) 



(T + R)" = Env^ r ^" ( 42 ) 

n=0 
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Thus, the average and variance of the occupation numbers are 

(on) = ^ = T (45) 

(Sol) = (n V n)2 = T(1 " T) (46) 
F(f) = 1-T = F if (47) 

The subscript "if" stands for the independent Fermionic case. The Fano 
factor is less than unity, which is smaller than the Poissonian case. Hence, 
sometimes this is called subpoissonian shot noise. Just to show the full 
expression, the shot noise of the independent Fermionic wave packets is 

S(f) = 2q{o n )F (I i )G(f) = 2qT(l-T)G(f)(I i ) (48) 

Let us denote the shot noise for the independent fermionic case as S{f(f). 

S tf (f) = 2qTF if G(f)(I i ) (49) 

2.3.3 Decomposition of Fano factor 

Let us decompose the Fano factor to make its structure clearer. 

t? ( £\ - ( 6 °D , (SOn+mSOn) , s. 

[°n) m> > \O n ) 

(Sol) . (So n+m SOn) , , 

= -. — r + / ; — ; cos(mujT) 

(o n ) ^ (On) 
rn>0 

($°l) fi , (So n+m 5o n ) 



0On) {, . {dOn+mdOn) , s 



(°n) L f^n ( S ° 
n+n%S 

TsoJ) 



m>0 



, (So n + m So n ) 
1 + cos (™^) 



m>0 

= F A(f) (/>0) (50) 

The first term describes the correlation within an event and the second term 
describes the correlation among events. Let us simply say that the first and 
second term corresponds to the intra-event and inter-event correlation, re- 
spectively. While the intra-event correlation gives a whilte spectrum for the 
Fano factor, the inter-event correlation induces the frequency dependence. 
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2.3.4 Low frequency shot noise and effective charge 

At low frequencies / << r _1 , it is clear that the low- frequency component 
of the shot noise is insenstive to the details of the wave packet shape as 
long as the J dtg(t) = 1, namely, its DC component is unity. Furthermore, 
cos(mojT) that comes into the Fano factor goes to 1 in the low frequency 
limit. In short, cos{mojT) — > 1 and G{f) — > 1 in the limit of / — > 0. 

A(/-X>) = [l+X; <%sfe> ] - Ao (51) 
L „^ (So 2 n ) J 

Note that Ao = 1 in the independent Fermionic case. Therefore, the low- 
frequency shot noise is as follows. 

S(f^0) = 2qTF o A (I i ) (52) 

Now we are at the stage to introduce the concept of effective charge. The 
effective charge is a covenient picture to capture the tunneling process intu- 
itively instead of talking about the exact values of Fano factor. It interpretes 
the tunneling process such that Fermionic wave packets with charge q* in- 
dependently tunnels across the partitioner. Namely, 

S(/->0) = 2q*TF if (Ii) (53) 

This is the definion of the effective charge. We can rewrite it as follows. 

S( f -> 0) F A a 

5 = sm^W)^ - (54) 

As mentioned before, the effective charge is effective since the details of wave 
packet shape is irrelevant in the low- frequency limit. There are two ways 
to modify the effective charge q* from q, which are intra- and inter-events 
correlations. To look at the effect of Fq and Ao separately, we will consider 
two cases, Fq = Fif and Aq = 1 in details. 
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2.3.5 Effective charge via correlation among tunnelings 

In the independent Fermionic case, the charge of the tunneling entity cannot 
be deviate from q. We assume here, for simplicity, that Fq = Fjj, which 
corresponds to systems with mere electrons as only available excitations (no 
exotic fractional quasiparticles). Thus, the effective charge is now 

q* = Aoq (55) 

Let us take several examples to capture the idea of inter-event correlation. 
Note that the present method does not provide the microscopic mechanism 
how such correlation is induced. It rather allows to compute the effective 
charge for a given set of correlations. Hence, one can also come up with a 
list of possible correlations to yield the observed effective charge. To resolve 
among the possible types of inter-event correlations, one has to measure 
the full power spectrum. Suppose there is a correlation in only neighboring 
events and no further correlation. 

(So n+1 So n ) = ^{5o 2 n ) (56) 

-> Ao = 1 + f = | (57) 

There are two main processes where one wave packet transmits alone and 
two consecutive wave packets transmit together. The two-particle process 
looks as if they are bunching. Sometimes it is called particle bunching. The 
effective charge is §<?. Next, suppose there is an anti-correlation in only 
neighboring events. 

(5o n+1 5o n ) = ~(5o 2 n ) (58) 

-+ A = 1-| = i (59) 

Thus, the effective charge is ^q. Let us look at another example with a 
further inter-event correlation. 

(5o n+2 5o n ) = ^(5o 2 n ) (60) 
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(5o n+1 5o n ) = ^(5ol) (61) 

-> A = l + ^ + ^ = 2 (62) 

The effective charge is 2q. The last example also gives rise to the effective 
charge of 2q but with different inter-event correlations. 

(5o n+3 5o n ) = ^(5ol) (63) 
(5o n+2 So n ) = I (So 2 n ) (64) 



(8o n+1 So n ) = \{5o 2 n ) (65) 



- Ao = ^+\+\+\ = 2 (66) 

In actual experiments, one often faces irrational effective charges even in 
systems with, presumably, mere electrons. Quantum dots are good example. 
With multiple states within the bias window or cotunneling by higher order 
processes can enhance the effective charge from mere q. In order to reveal 
whether the effective charge is the elementary charge of the system or not, 
one should take T ~ or T ~ 1 where the single particle process dominates 
and the distribution reduces to the Poissonian case. 

2.3.6 Effective charge via uncorrelated multiple occupation 

Suppose now each slot of an ordered train accomodates up to K wave packets 
(A; = 0, 1, 2, K) and there is no inter-event correlation (Ao = 1). 

<f = (67) 

Let us denote the probability that k wave packets transmit together as 
with X^fc=o = 1- The distribution that we employ now is the multinomial 
distribution. There are NK wave packets over the train and tunneling 
events that k wave packets transmit together. Hence, Yl,k=o n k = N. 

K 



fc=0 . 



-I N 

Y , f -T"°T^---T^ K (68) 

(no,ni,—,n K ) 
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Given the distribution, the expectation values of n k and n\ are 



<n fc > = n^L 



(4) = n 



K 

E r * 

s=0 

d 



—T - 



N 



K 



= NT k 

N 



(69) 



s=0 



= (NT k ) 2 + NT k (l-T k ) (70) 



The expectation values of the mean and variance of the occupation o n and 
the Fano factor read 

K K 

X 



(On) 

(So 2 n ) 



j K K 
— J2k{n k ) = J2kT k 



r=0 
K 



k=0 



^f> 2 «> = Efe 2 T fc (l-T fc ) 



fc=o 

Ef =0 fc 2 r fc (i-T fc ) 



k=0 



(71) 
(72) 
(73) 



l>2k=o kTk 

One may define the normalized transmission probability T nrm and the nor- 
malized Fano factor as 

K 



F, 



O.nrm 



(On) 

K 
K 



= E 



K~ 



k=0 

Ek=o(Tt) 2 Tk(l-Tk) 



(74) 
(75) 



EK k_rp 
fc=0 K ±k 

Tnrm is nothing but the transmission that is measured in actual conductance 
measurements. Note that the initial current is = — 2 and the number of 



events per unit time is 



= 2q(o n )F ^G(f) = 2KqT nrm F 0ynrm (I l )G(f) (76) 

Alhough the multiple occupation case sounds unphsycial, the edge state of 
fractional qunatum Hall effect, e.g. the filling factor 1/3, may be modeled 
as the case of K = 3 with q = |. Let us consider the low- frequency part 
of the shot noise given in the last section. S(f ps 0) = 2KqT nrm Fo tJirm (Ii) . 
The effective charge can be defined with the shot noise of the indepen- 
dent Fermionic wave packets with transmission! probability T nrm , which is 
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2qT nrm (\ T nrm )(Ii) . 

* F KF 0nrm 
Q = p 1 = -B 1 ( 77 ) 

r if,nrm r if,nrm 

, where F iftnrm = 1 - T nrm and F = KF 0tnrm . 

2.3.7 Effective charge via correlated multiple occupation 

Though we do not go into details, for the completion, we show the expression 
of the effective charge in the case of multiple occupation with correlation 
among the tunnelings. 

* KFo nrm Ao 

q = ~T~ f q 

2.3.8 Case of uncorrelated double occupation 



(78) 



Here, we consider the case of K = 2 explicitly. 

(1) : (T ,r 1 ,T 2 ) = (i?,r,0) 

((o n ),(5o 2 n ),T nrm ,F ) = (T, T(l-T), f, 1 — T) 

q*(T) = \^q, q*(T^0) = q (79) 

(2) : (T ,T l ,T 2 ) = (R,^,^) 

((o n ),(5o 2 n ),T nrm ,F ) = (3f, 5f(l-f), |f, |(1-|)) 

Q*(T) = q*(T^0) = \q (80) 

(3) : (To,Ti,r 2 ) = (i2,0,T) 

((o n ),(5o 2 n ),T nrm ,F ) = (2T, 4T(1-T), T, 2(1 - T)) 

g*(T) = 2q, q*(T^0) = 2q (81) 

(4) : (T ,T 1 ,T 2 ) = (R,aT,(l-a)T), (0 < a < 1) 
((on), (H 2 t )) = ((2 - a)T, (4 - 3a)T - (4 - 8a + 5a 2 )T 2 ) 
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Here, we look at them a bit more carefully. 



T nrm = -(2-o)T (82) 



4 -3a/ 4-8a + 5a 2 T 
2-fll 4 - 3a 



) (83) 



A Q 1 _ 4-8a+5a z m 

9(T ' fl) = ^• l-l(2-a)T - g (84) 
g*(T^0,a) -> (85) 

It is interesting to notice that a = | yields q* = ^(g + 2g) = |g. The equal 
mixing of the one-particle and two-particle tunneling events (a = ^) does 
not simply leads to the effective charge right at the mid value q* = |g. This 
is because the variance goes as quadratic in the tunneling charge. 



2.3.9 Case of uncorrelated triple occupation 

Here, we consider the case of K = 3 explicitly. (1): (To, T\, T2, T3) = 
(R, aT, IT, {I- a- b)T) 

((o n ),(5o 2 n ),T nrm ,F ) = (T, T(l-T), ±T, (1 - T)) 

T nrm = -(3-2a-6)T (86) 
(9 - 8a - 56) - (9 + 10a 2 + 136 2 - 9a - 96 - 9ao)g ^ 



3 - 2a - b 

, _ (9+10a 2 +13b 2 -9a-9fc-9ab) ^ 

* (I> ' 6) = 3-2,-i, 1 - 1 (3 - 2a - b)T " (88) 

2.3.10 Case of multiple independent Fermionic conductors 

Suppose now there are multiple (ID) conductors participating in the trans- 
port. However, to keep it simple, let us assume they are just independent 
from each other (no wave packet swapping) . Let us index the conductors as 
I = 1, 2, ■ ■ ■ , L. The full expression is quite straightforward. If we suppose 
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all the wave packet shapes in the conductors are identical, Gi(f) = Gf. 

L 

S(f) = X>( On) ,)F 0) , <£><?,(/) (90) 
1=1 

L 

= £2^,(1 -T,)<ii>G,(/) (91) 
l=i 

L 

= [£2<ffi(l-T,)](Ii)G(/) (g2) 
i=i 

One can think of a model system as a quantum dot or single molecule bridge 
where multiple channels for the transmission can be available. For that a 
more careful consideration such as assigning different Gi(f) and (In). The 
index I is for each transport channel. In the Landauer's paper, he considered 
the shot noise as a function of the DC bias eV. First, he split this energy 
window into many tiny energy segments AE, calculated the shot noise, and 
then summed over the entire energy with occupations follwoing the Fermi- 
Dirac distribution. Regarding each segment as a conductor, we can recover 
and even extend the Landauer's approach with the arguments given so far. 

2.3.11 Case of fluctuating tunnel barrier 

We again consider the independent Fermionic case. However, now the tun- 
neling probability is fluctuating in time. We start with the simplest case. 
Namely, it takes two transmission probabilities T\ for < t < N\t and T 2 
for Nit <t<{Ni + N 2 )t, where N = iVi + N 2 . Important to note that, at 
the end, Ni and N 2 should also go to infinity in taking the limit of 2? — > 00 
with keeping the ratio Ni/N and N 2 /N. The distribution of the transmitted 
wave packets is 

( /Vi 

711=0 



N 2 



I" ^ rpn-z nAT 2 -n 2 

[L, (N 2 -n 2 )\n 2 \ l2 K " 

712=0 
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Thus, the number of the transmitted wave packets n is n = rii + n 2 . We 
would like to calculate the expectation value and the variance of n. First, 
the expectation value of n is 



(n) = (ni) + (n 2 ) = + N 2 T 2 (94) 

(On) = -Jf = J^T 1 + —T 2 (95) 

THen, to calculate the variance, the expectation value of n 2 is 

(n 2 ) = {n\) + 2(mn 2 ) + (n 2 > 
= (n 2 )+2(m)(n 2 ) + (n2) 
(,5n 2 ) = ((5n 2 ) + (Jn 2 ,) (96) 

The decomposition of (nin 2 ) into (ni)(n 2 ) is simply because the binomial 
distribution assumes the independent events. 

= ^ = §T l( l-T l) + §T 2 (l-T 2 ) (97) 

By induction, the extension to more fluctuating case is trivial. For the 
case of M tunneling probabilities, the distribution function, the expectation 
value, and the variance is 

M M AT m j 

nor™ = n [ e (N _Zvn ^ mR ^ nm ] ^ 

m=l m=l n m =0^ Vm n "^ n ™- 

1 M 

(On) = JjY, NmTm (") 
m=l 

1 M 

= (ioo) 

m=l 

Now, we can take two ways in taking 2? — > oo. One is to keep he ratio N m /N 
and the other is to take M to the infinity. For the former, the result given 
above is enough. However, for the latter, we need to know the distribution 
of T m to obtain the proper expectation values. 
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3 Practical Implementation 



3.1 Finite Bandwidth Detector 
3.1.1 Center frequency = DC 

In practice, the detector has a finite time resolution, At, with a specific shape 
of time-window, Z(t). Z(t) monitors the current reaching the detector, I(t) 
and smears the original singal. The singal after the detector, Init), becomes 
the convolution of the original signal and the time window, I(t) * Z(t). In 
the Fourier space, the detector acts as a filter with a frequency response 
Z(f) = &\Z(t)\ centered at / = 0Hz and a bandwidth A/. 

I D {t) = I(t)*Z(t) (101) 
&[Z{t)\ = Z(f) (102) 

In most electronic systems, At m (A/) -1 and At >> r. Of course, the 
exact relation depends on the detail shape of the filter. 

Given the realistic detector, what is the measured shot noise power spec- 
trum, S D (f)7 

S D (f) = J?[C ID (t)] + c.c. (103) 

C D (t) = (SI D (t')6I D (t' +t)) (104) 

5I D (t') = J dt"Z(t")5I(t' -t") (105) 

SI(t'-t") = Y,o k g(t'-t"-t k ) (106) 

k 

Using above fomula, let us calculate the correlation function and the power 
spectrum explicitly. 

C D (t) = lim -1 I dt'6I D {t')5I D (t' + t)) 

= J im |f Yl ° k ° k ' I dt ' I dt " Z «)9(t' - t'{ - t fe ) I dt'iZ{dt'i)g{t + t' - 4 - t fe 
: ~ k,k' ^ J J 

= J im §F S ° k ° k ' I dt ' I dt'lZ^'lW - t'{ - t k ) I d%Z(4)g(t + 1' - 4 - t k ,) 
' J k,k' J J J 
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q 2 * f 
« — y~]((5o n+m (5o n ) dt'u(t')u(t + t f + mr) 

m 

= g(/ J )^(fc„+ m 5o n )C/(t + mr) (107) 

m 

Here, we have introduced two functions: 

u(t) = y dt'Z(t')g(t - t') (108) 
I7(t) = y dt'u(t')u(t' + 1) (109) 

Their Fourier transforms are 

J?Ht)] = u(f) = Z(f)~g(f) (110) 
W*)] = U(f) = u(f)u*(f) = \Z(f)\ 2 \~g(f)\ 2 (HI) 

Therefore, the measured shot noise Smif) is given by 

S D (f) = 2q(I i )J2(S°n+mSo n )cos(muT)U(f) (112) 

m>0 

= 2 gJ F|Z(/)| 2 [/(/)( 0n )(/,) (113) 



= \Z(f)\ 2 S(f) (114) 
3.1.2 Center frequency > DC 

Suppose that we shift the center frequency to fa. The consequence is trivial 
in the frequency domain. 

S D (f) = \Z(f-fa)\ 2 S(f) (115) 

The effect of the frequency shift in the time domain is also simple. 

3?[Z(f-fa)} = e^^Z(t) (116) 



Re 



&[Z(f-fa)\ = cos(2irfat)Z(t) (117) 



Thus, the effect of shifting the center frequency on the detector function 
Z(t) is merely modulating it with cos(2ir fat). 
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3.2 Cross Correlation 

3.2.1 General form of auto correlation 

The shot noise formula given above is as follows. 

C(t) = {5I(t')5I(t' + t)) (118) 
S(f) = 2\*[C{t)]\ (119) 

This is nothing but the auto correlation function of the current fluctuation. 
We name this current as current 1, Ii(t) = qigi(t) * h\(t) with tik = kr\. 
Then, we may denote it as follows. 

C(t) C ll {t) = (5h{t')5h{t' + <)) (120) 
S(f) S n (f) = 2\^[C u (t)]\ (121) 

Let us also define the following. 

G(t) G n (t) = J dt' gi (t') 9l (t + t') (122) 

(Oln) 



F =► Fn(f) = ^ {SOln . +mS ° ln) cosimum) (123) 



m>0 

One can also include the detector as previously discussed. 

U(t) Un(t) = J dt'ui(t')ui(t + t') (124) 

ChHt) = (SI m (t')SI D1 (t' + t)) (125) 
Sh\f) = = 2\*[C#{t)]\ (126) 

So far, this is nothing but a change of notation. The benefit from this 
generalization becomes clear in the followings. 

3.2.2 General form of cross correlation 

We may generalize this formula by considering the cross correlation func- 
tion between two kinds of currents. Namely, I\(t) = q\g\(t) * h\{t) and 
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h(t) = q292(t) * h 2 (t), where hi(t) = Ylk=i°^(t - t lk ) and h 2 (t) = 
°2fc'<5(* - hk')- As before, t lk = fen and t 2k > = k'r 2 . q\ / q 2 means 
two currents consists of two different charged entities, n / t 2 means the 
two currents consists of slots with different duration. However, they should 
satisfy N\T\ = N 2 t 2 . To be more concrete, two currents are coming from 
different contacts which are voltage-biased differently. o\k can obviously be 
different o 2 k> in general. The cross correlation of the current fluctuations of 
these two currents is defined as follow. 

Ch 2 (t) = (SI D i(t')6I D2 (t' + t)) (127) 

The detailed calculation follows: 



1 r J i z 

C%(t) = lim — / dt'5I D1 {t')5I D2 {t' + t) 

iVi N 2 

= ^fim^ ^2 ^2 ^°ifc^°2fc' / dt'ui(t')u 2 (t + t' + t lk -t 2 k>) 



k=i k'=i 

Ni N 2 

= Y\m^ —gT y $Qik8o 2k 'Ui 2 (t + ti k ~ t 2k >) (128) 
fc=ifc'=i 

, where, obviously, Ui 2 (t) = f dt'ui(t')u 2 {t + t'). With this maximum gen- 
erality, it is difficult to proceed much further. We have to take some simpli- 
fications. However, in general, one is not interested in cross correlating two 
signals prepared completely differently. Let us say that our main playground 
is the mesoscopic electric transport, one is oftern interested in cross corre- 
lating two currents after partinining a current or two currents with same 
voltage-biases such as in collision experiments. Hence, without losing much 
of generality, it can be simplified once we take the case n = t 2 = r, thus, 
Ni = N 2 = N. 

N 
k,k 

N-l ra=[|m|+l:JV],m<0 



Cd(*) = i im y Soi k So 2k 'Ui 2 (t + t k - t k ' 

k,k'=l 



lim ^2 ^2 Soi n+m 5o 2n Ui 2 (t + mr) 

m=-(N-l) n=[l:JV-m],m>0 
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mr) 



~ "^y^(8o ln+m 5o2n}Ui 2 (t + rnT) 

T *■ — ' 

m 

Q1Q2 n rj r-A (5oin+m5o 2 n} TT , 

, Iqi(oin) Iq2(02n) (So ln+m So 2n ) TJ ,. , \ 

= ^^K¥E ( 5 °^rn60 2n ) ^ + (l2g) 
m V(°ln)(°2n) 

Its power spectrum follows 

Sh 2 (t) = 2^^! ^±f^ c ,(ma;r)£/ 12 (/) (130) 

m > V(°ln)(°2n) 

= 2y/q^y/{h){I 2 )F 12 (f)U 12 (f) (131) 
If q 1= q 2 = q, = (I 2 ). 

* 

C%(t) = q(Ii)J2( S °in+mSo 2n )U 12 (t + mT) (132) 
m 

* 

S^t) = 2g(ij) ^ (<5oi n+m (5o2n)cog(majr)i/i2(/) (133) 



m>0 



3.2.3 Cross correlation of electrons from single splitter 

Let us make an exercises on this cross correlation. The simplest case is 
cross correlation between the transmitted and the reflected currents from a 
quantum point contact (QPC) in the filling factor 1 of quantum Hall system. 
A noiseless current Ii(t) = (ij) is impinged on the QPC. q\ = q 2 = e and 
wave packets obey the Fermi statistics. The occupation of each slot is, hence, 

02k = 1 - oik => So 2k = -So lk = -So k . (134) 



Therefore the cross correlation of two currents is 

CbHt) = -e^WMjZ {6 °; +m *° n) U 12 (t + mT) (135) 

= -eVm^t {6 ;; +m ' 0n) Un(t + rnr) (136) 

m V(°ln>(02n> 

= -Ch\t) (137) 
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The power spectrum of this cross correlation is 




cos(muT)U 12 (f) (138) 



(139) 



= -2ey / {oin)(o2n)F 12 (f)U 12 (f)(I i ) 



(140) 



-2e^T(l-T)F 12 (f)U 12 (f)(I l ) 



(141) 



- -2ey/T(l-T)Fn(f)Uu(f)(Ii) 



(142) 



= -sh\t) 



(143) 



Here, we used {o\ n ) = T, (o 2n ) = 1 — T. Hence, making a cross correlation 
on the transmitted and reflected currents gives the same result as the auto 
correlation other than the sign. However, in reality, noises of the two de- 
tectors are uncorrelated. By averaging the uncorrelated noises, one can get 
a better signal to noise in the measurement of the shot noise compared to 
detecting it with one detector given the same averaging time. 

3.2.4 Case of independent tunneling 

For the case of non-correlated tunneling events, (8o n+m 8o n ) is finite only for 
m = 0. Remember (50%) = T(l — T) for independent tunnelings. 



Remember that G± 2 (t) = Gn(t) = G 22 (t) since the two currents are origi- 
nating from the same current. Therefore, the power spectrum of the cross 
correlation reads 



, which is nothing but the shot noise of electron currents with negative sign. 
This negative sign originates from the conservation of current and does not 
prove anti-correlation among Fermions. 



C v \t) 



eTil-T^Uuit) 



(144) 



S 12 (t) = -2eT(l-T){Ii)U u (f) = S n (t) 



(145) 
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3.2.5 Effect of time-delay on cross correlation 

In reality, the distances between the tunneling barrier and each detector can 
differ. This means that we cross correlate signals at different times. Let 
us name it as time-delayed cross correlation. This time delay can be also 
intentionally introduced using a proper electronic component. We can easily 
implement this time delay t d < as follows. 

C%{t + t d ) = (6I 1 (t')6I 2 (t' + t + t d )) (146) 

The correlation function and the power spectrum are easily computed. 

ChHt + t d ) = V^Vih){h) E ^ 0l "+- fo2 "> Ul2 (t + t d + mr) (147) 

m>0 V \°ln/V \°2n/ 



= 2^^^/{h){h)F 12 (f,t d )U 12 (f) (148) 

= Wi sl2(m (149) 

It is clearly seen that the time delay t d does have an effect on the cross 
correlation by shifting the phase. To understand it better, let us again take 
the case of independent tunnelings, namely, keeping only m = in the sum. 

(5oi n 5o 2n ) 



S%(f,t d ) = 2^^!^) ^11 p= cos(ut d )U 12 (f) 

V (°ln)V \02n) 

= -2eT(l - T)U 12 {f)co8(wt d ){Ii) 

= cos(ut d )S%{f,0) (150) 

For the center frequency / =lMHz, suppose there is a 3-meter difference 
in the cables for the detectors. Disregarding the dielectric constant of the 
coaxial cable, t d = 3xl Q? m/ / s = lOnsec. Then, cos(uot d ) ps 0.998. (Of course, 
in reality, the coaxial cable adds extra stray capacitance.) Or, suppose 
that the two cables are equal in length and that the path-length difference 
between the tunnel barrier and the detectors is 1mm. For the Fermi velocity 
of 10 5 m/s, t d = = lOnsec and cos{ut d ) w 0.998. 
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